Search for Cosmic-Ray Antideuterons 
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We performed a search for cosmic-ray antideuterons using data collected during four BESS balloon 
flights from 1997 to 2000. No candidate was found. We derived, for the first time, an upper limit of 

I. 9xl0 -4 (m 2 s sr GeV/nucleon) -1 for the differential flux of cosmic-ray antideuterons, at the 95% 
confidence level, between 0.17 and 1.15 GeV/nucleon at the top of the atmosphere. 

PACS numbers: 98.70.Sa, 95.85.Ry, 96.40.De, 97.60.Lf 
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The possible presence of various species of antimatter 
in the cosmic radiation can provide evidence of sources 
and processes important for both astrophysics and ele- 
mentary particle physics. 

For example, discovery of a single antihelium nucleus 
in the cosmic radiation would offer clear evidence for a 
baryon symmetric cosmology. Despite extensive and on- 
going searches, none has ever been found 0,0- 

Similarly, the spectral form and magnitude of the an- 
tiproton (p) spectrum could provide evidence for a num- 
ber of possible primary sources, including evaporating 
primordial black holes (PBHs) 0, Q and annihilating 
neutralino dark matter J^, 0] as well as a baryon symmet- 
ric cosmology. However, recent results from the BESS ex- 
periment UH;I3 imply that most of the p's in the cosmic 
radiation are not from primary sources, but rather are 
secondary products of the energetic collisions of Galactic 
cosmic rays with the interstellar medium. The data do 
suggest that below ~1 GeV there is an excess of p's above 
expectation from a purely secondary origin, but the situa- 
tion is far from clear. Model calculations of the secondary 
spectrum still have ambiguities 0, 01 El El H3| and 
statistical errors of the currently measured low-energy 
p spectrum are not small enough to provide clarity. The 
accuracy of both calculations and measurements needs 
substantial improvement. 

While antideuterons (d's) have never been detected 
in the cosmic radiation, they can be produced by the 
same sources as p's and may be of both secondary or 
primary origin, with the latter providing evidence for 
sources such as PBHs and annihilating neutralino dark 
matter. The low energy range below ~1 GeV/nucleon 
offers a unique window in the search for cosmic-ray pri- 
mary d's because it has a grea tly reduced background 
from secondary d's EI Ilfij. as compared with sec- 
ondary p's. Thus, the unambiguous detection of a single 



d below ~1 GeV/n would strongly suggest the existence 
of novel primary origins. Hence, cosmic-ray d's have an 
advantage over cosmic-ray p's as a probe to search for 
primary origins. 

In this paper we report on a search for d's carried out 
with four balloon flights of the BESS instrument from 
1997 to 2000. Using data from these flights, we report 
for the first time an upper limit on the differential flux 
of cosmic-ray d's and discuss this result in the context of 
expectation from evaporating PBHs. 

The BESS detector was designed El El and devel- 
oped 0| as a high-resolution spectrometer with the large 
geometrical acceptance and strong particle-identification 
capability required for antimatter searches. A uniform 
magnetic field of 1 Tesla is generated by a thin su- 
perconducting solenoid. The field region is filled with 
tracking detectors consisting of a jet- type drift cham- 
ber (JET) and two inner drift chambers (IDCs). Track- 
ing is performed by fitting up to 28 hit points in these 
drift chambers, resulting in a rigidity (R) resolution of 
0.5% at 1 GV. The upper and lower time-of-flight scin- 
tillator hodoscopes (TOFs) measure the velocity (j3) and 
the energy loss (dE/dx). The time resolution of each 
counter is 55 ps, which yields a 1/(3 resolution of 1.4%. 
A threshold-type Cerenkov counter with a silica-aerogel 
radiator (n=1.03 in 1997 and n=1.02 thereafter) can re- 
ject e~ / '[iT events, which can be backgrounds for the 
detection of p's and d's, by a factor of ~10 3 . 

Four balloon flights were carried out in northern 
Canada, 1997 through 2000: from Lynn Lake to Peace 
River where the geomagnetic cutoff rigidity ranges from 
0.3 to 0.5 GV. Data for the d search were taken for 
live time of 15.8, 16.8, 27.4, and 28.7 hours in 1997, 
1998, 1999, and 2000 respectively, at altitudes around 
36 km, corresponding to ~5 g/cm 2 in residual atmo- 
spheric depth. The data acquisition sequence was ini- 
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tiated by a first-level trigger, which is generated by a 
coincidence between hits of the top and bottom TOFs 
with a threshold set at 1/3 of the pulse height from min- 
imum ionising particles. In addition to biased trigger 
modes 0,^1 enriching negatively charged particles, one 
of every 60 (30 in 2000) first-level triggered events were 
recorded as unbiased samples. 

The concept of the off-line analysis is similar to that 
used for the p selection described in Ref @- For events 
of both negative and positive curvature, same selections 
were applied to detect clear d and deuteron (d) candi- 
dates. The selected d's were used to estimate selection ef- 
ficiencies for d's. At the first step, we selected events with 
a single downward-going, passing-through track which is 
fully contained inside the fiducial volume with restricted 
number of TOF hits, in order to reject interacted events 
as well as albedo particles. At the second step, in order 
to eliminate backgrounds such as large-angle scattered 
events by ensuring good quality of R and /3 measure- 
ments, we applied several cuts on tracking and timing 
measurement quality parameters such as: (i) the number 
of used hits and the reduced x 2 of the trajectory fitting, 
and (ii) the consistency between the JET track, hits in 
the IDCs, and the TOF timing information. 

In order to identify d's, dE/dx measurements inside 
the TOFs and the JET were required to be consistent 
with d's as a function of R. In addition, the Cerenkov 
veto was applied to reduce the e~/p~ background con- 
tamination. Thereafter, the mass of the incident particle 
was reconstructed using the measured j3 and R. Figure^] 
shows the 1//3 vs. R plots of the events which survived all 
the above selections (for 1997 - 1999, only the negative 
rigidity events are shown). The d selection region was 
determined by the mirror position of the d band, which 
was defined to have a uniform selection efficiency of 99%. 
In order to avoid the contamination (or misidentification) 
of p's, the region overlapped by the p band was excluded. 
The exclusion band was defined so that it had a uniform 
selection efficiency for p's and the possible p contamina- 
tion from the whole data set was just one event or less 
in an estimation using the 1 / (3 distributions of positively 
charged events. The contamination of other negatively 
charged particles (mainly e~ was estimated to be 
less than 0.1 event. The possibility of the spill over of 
positively charged particles into the negative side is neg- 
ligible in the considered rigidity region. In Fig. ^ (panel 
BESS00), candidates of protons (p's) and p's selected by 
the same procedure are superimposed. 

As shown in Fig.^ clear mass bands of p's, d's, tritium 
and p's can be seen. However, no d candidate exists 
within the expected selection region. 

Since no d candidate was found, we calculated the re- 
sultant upper limit on the d flux |ljj , which is given 
by: $ d - = N ohs /\SQ e tota] (1 - <5 sys )| min /T live /(£ 2 - E 1 ). 
The live time, Ti; ve , was directly measured by a 1 MHz- 
clock pulse generator and scalers throughout the flights. 
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FIG. 1: The surviving single-charge events in the data of 
each flight. The dotted curves define the d mass bands. The 
dotted vertical lines at ~-3.7 GV corresponds to E2 shown in 
Fig. El 



As the number of the observed d events, A^bs, we took 
3.09 for the calculation of the 95% C.L. upper limit. We 
did not consider the effect of the possible background 
contamination (<1 event), because the background es- 
timation still has an ambiguity since it is difficult to 
evaluate the amount of the tail distribution strictly. E\ 
and E 2 denote the energy range of the limit at the top 
of the atmosphere (TOA) . The energy measured at the 
instrument was traced back to the one at the TOA by 
correcting the ionisation energy loss. In order to ob- 
tain the most conservative limit, the minimum value of 
(Sfl etotai (1 — feys)) was used, where Sfl is the geomet- 
rical acceptance, etotai is the total detection efficiency, 
and S S ys is the total systematic uncertainty. The etotai 
can be written as e to tal = Etrig £1 £2 e P id £ aC c e a ir ■ The 
efficiency of the first step selection including the effects 
of inelastic interactions in the instrument (ei), the sur- 
vival probability in the residual atmosphere (£ a ir)i an d 
the SQ were calculated by the BESS Monte Carlo (MC) 
based on GEANT/GHEISHA. Since there are no exper- 
imental data of d interactions in material, we incorpo- 
rated the d in the code under the following assump- 
tions: (i) The inelastic cross sections of d can be esti- 
mated by scaling those of p usingan empirical model of 
hard spheres with overlaps I2H Eli , which is described 
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where a(Ai,A t ) is the cross section of an incident parti- 
cle with atomic weight Ai to a target with atomic weight 
At- (ii) When an inelastic interaction occurs, d is al- 
ways fragmented or annihilated, (iii) Other effects of 
energy loss, multiple scattering, bremsstrahlung and 8- 
rays are described as are those of other nuclei. This 
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FIG. 2: Effective exposure factors of each flight data and their 
sum. The systematic uncertainty is shown as a ±10% width 
band. E\ and E2 denote the energy range of the limit. 



hard sphere model is known to reproduce data on nu- 
clear interactions for various combinations of Aj and A t 
including light nuclei such as p/d 0, and antinu- 
clei including p/d (2^, 0- We adopted this model to 
estimate the a(d, A t ) scaling from the a(p, A t ) described 
in Ref. j^. The efficiency of the the second step selec- 
tion (£2) was estimated by using both the unbiased data 
and the BESS MC. The trigger efficiency (etrig) was ob- 
tained by using the unbiased data and a detector beam- 
test data |25j. The efficiency of particle identification 
(£pid) w as estimated using the unbiased d samples of each 
flight under the assumption that the d candidate should 
behave similarly to d except for deflection in the symmet- 
rical configuration of BESS. Typical values at 0.5 GeV/n 
are: £trig ~ 90%, e x ~ 60%, e 2 ~ 70%, £ pid ~ 98%, 
£air ~ 85%, and SQ ~ 0.25 m 2 sr. The probability of 
events without any hits or tracks by another accidental 
incident particle, £ aC c, was derived to be ~94% by sam- 
ples taken by the random trigger which was issued at 
once a second throughout the flights. 

FigureElshows the calculated effective exposure factor. 
The decrease of the factor at the low-energy side is mainly 
caused by the decrease of the geometrical acceptance, the 
decrease of mean free path through the detector, and the 
increase of large-angle scattering. The major reason for 
the decrease at the high-energy side is the decrease of £ p id 
due to the overlap of 1//3 distributions between d's and 
p's. The combined systematic uncertainty, which was es- 
timated to be <5 sys ~ 10% with less energy dependence, 
is also shown in the figure. Dominant systematic uncer- 
tainties were the uncertainties in the evaluation of £1, £2 
and £ a ; r , all of which were discussed using the MC sim- 
ulation. The energy range of E\ - E2 was chosen to be 
0.17 - 1.15 GeV/n, where the exposure factor is highest 
and has relatively little energy dependence. 

The resultant upper limit $ dT for 1997, 1998, 1999, 
2000, and the integrated flight data were calculated to 
be 9.8, 8.9, 6.9, 6.2, and 1.9 xl0~ 4 (m 2 s s^GeV/nT 1 
respectively (Fig. These are the most conservative 
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FIG. 3: Obtained upper limit on the d flux in comparison 
with the PBH-d spectrum calculated in this paper. A p spec- 
trum measured by BESS and some theoretically predicted 
spectra of p and d are shown for reference. 



limits with no assumptions on the d spectrum shape. 
If we assume a uniform d energy spectrum and use a 
mean inverse exposure factor, the summed upper limit 
was evaluated to be 1.6xl0~ 4 in the same energy range, 
and 1.4xl0~ 4 in the range 0.13 - 1.44 GeV/n where the 
upper limit is minimized under this assumption. Since 
our detection efficiency of d is less dependent on the en- 
ergy, the upper limit is less dependent on the assumption 
of the d spectrum shape. In the following discussions, we 
use the most conservative one (1.9xl0 -4 ). 

As described in Ref. Q, only PBHs that are close to 
explosion and exist within a few kpc of the solar sys- 
tem can contribute to the observed flux. Therefore, the 
d upper limit leads directly to the upper limit on the 
explosion rate of local PBHs, 7?.pbh- In order to ob- 
tain 7?.pBH from the d upper limit, we have calculated an 
expected PBH-d spectrum through the following steps: 
(i) the emission rate of particles from PBHs, (ii) the 
fragmentation rate to form d's, (iii) the source spec- 
trum, (iv) the propagation process, and (v) the effect 
of the solar modulation. The calculations except for (ii) 
were based on the calculation of PBH-p spectrum de- 
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scribed in Ref. |3j. The step (ii) was performed by using 
the frequently-used "coalescence model" (e.g. Ref. [26jk 
According to this model, the production probability of 
d'a in momentum space, d 3 rtj/dp 3 , can be expressed as 
the product of those of p's and antineutrons: 7 
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) ~ t^Po (7 li^) 2 : wh ere p is 



the "coalescence momentum" which must be determined 
from experiments. We assumed po = 130 MeV/c from 
the data of d production in e + /e~ annihilation |27|. The 
solar modulation in step (v) was estimated by using the 
numerical solution of the spherically symmetric model 
proposed by Fisk [2^]. The solar modulation parame- 
ter, (j), was determined to fit the p spectrum measured in 
the same BESS flights as: 500, 610, 648, and 1334 MV 
in 1997, 1998, 1999, and 2000 respectively @, @. The 
calculated PBH-d flux for 1997 on the assumption of 
ftpBH = 2.2 x 10" 3 pc^yr^ 1 is shown in Fig-EJ The dif- 
ference between our PBH-d flux and the one in Ref. |2^| 
mainly comes from the different assumptions on the prop- 
agation model, similar to the case of the difference in the 
PBH-p spectrum between Ref. ; 3j and Ref. 4] that is 
described in Ref. 0. 

Here, we place the upper limit on the 7?.pbh to be 
1.8 x 10° pc _3 yr _1 , which is five orders of magnitude 
more stringent than the sensitivity for 50-TeV 7-ray 
bursts [33- The limit on TIpbh leads to an upper limit 
on the density parameter of PBHs in the Universe, f2pBH> 
to be 1.2xl0 -6 . The initial mass spectrum of PBHs was 
assumed to have a — | power-law form, and the PBH 
spatial distribution was assumed to be proportional to 
the mass density distribution of dark matter within the 
galactic halo Q. 

As a conclusion, we have searched for cosmic-ray 
d's with the BESS flight data obtained between 1997 
and 2000. No d candidate has been detected. Wc 
placed, for the first time, an upper limit of 1.9 x 
10" 4 (m 2 s sr GcV/n)-! (95% C.L.) on the differen- 
tial flux of cosmic-ray d's in an energy range of 0.17 
- 1.15 GeV/n at the top of the atmosphere. In conse- 
quence, we derived an upper limit of 1.8 x 10° pc~ 3 yr _1 
(95% C.L.) on the explosion rate of local PBHs and an 
upper limit of 1.2 x 10~ 6 (95% C.L.) on the density pa- 
rameter of PBHs. 

These upper limits regarding PBHs are two orders of 
magnitude looser than those derived from the p flux Q . 
However, further sensitive searches could push down the 
limits from gTs below the ones from p's, because the low- 
energy range has a greatly reduced background from sec- 
ondary d's. Astrophysical consequences of our d search 
will motivate further sensitive searches for d's as well as 
further advances in the physics of primary origins, in con- 
nection with cosmology and elementary particle physics. 
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